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Abstract 

In this article, we study the light-flavor scalar and axial-vector diquark states in 
the vacuum and in the nuclear matter using the QCD sum rules in an systematic way, 
and make reasonable predictions for their masses in the vacuum and in the nuclear 
matter. 
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1 Introduction 

We can rephrase the scattering amplitude of one-gluon exchange into an antisymmetric 
antitriplet 3c and an symmetric sextet 6c in the color-space, 

A" \ / A" \ 1 1 

yI (y) ^ --^{^jAi - ^ik^ji) + -^i^jAi + ^ik^ji) , (1) 

where the A*^ is the Gell-Mann matrix element, and the i,j and k, I are the color indexes of 
the incoming and outgoing quarks respectively. The attractive interaction in the antisym- 
metric antitriplet favors the formation of the diquark states in the color antitriplet, while 
the most stable diquark states maybe exist in the color antitriplet 3c, flavor antitriplet 3f 
and spin singlet 1^ channels due to Fermi-Dirac statistics [1]. We can take the diquarks 
as basic constituents to obtain a new spectroscopy for the mesons and baryons [21 [3], and 
the diquark states play an important role in many phenomenological analysis [HIS]. For 
example, we usually take the nonet scalar mesons below 1 GeV as the tetraquark states 
consist of the scalar diquark states [qq]-^^ and [qqjs^ in the relative S'-wave [5], and study 
the octet and decuplet baryons as the quark-diquark bound states [6j. 

The QCD sum rules is a powerful theoretical tool in studying both the in-vacuum and 
in- medium hadronic properties |7j, and has been applied extensively to study the properties 
of the in-vacuum hadrons and the in-medium light-flavor hadrons and charmonium states 
[SlIHlllQ]. In the limit niu = md 0, the in-medium nucleon mass can be related with 
the in-medium quark condensate {qq)pj^ through the simple relation = —YT{qq)pN^ 
where the is the Borel parameter. It is interesting to study the diquark states in the 
nuclear matter, as they are basic constituents of the baryons and play an important role 
in understanding the strong interactions and the relativistic heavy ion collisions. The 
in-medium baryon properties will be studies by the CBM (compressed baryonic matter) 
and PANDA collaborations unto the charm sector at the upcoming FAIR (facility for 
antiproton and ion research) project at GSI (heavy ion research lab) |1H I12j . The in- 
vacuum light and heavy scalar and axial-vector diquark states have been studies with the 
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QCD sum rules \13\ [T^ [15], we extend the previous works to study the in-medium mass 
modifications of the light-flavor diquark states. 

The article is arranged as follows: we derive the QCD sum rules for the light-flavor 
scalar and axial-vector diquark states in the vacuum and in the nuclear matter in Sect. 2; 
in Sect. 3, we present the numerical results and discussions; and Sect. 4 is reserved for our 
conclusions. 

2 The scalar and axial- vector diquark states with QCD 
Sum Rules 

We write down the two-point correlation functions Il{p) and H^uip) in the nuclear matter, 

n(g) = i I d'xe"^^^o\T{J{x)jHo)}\^o) , 
n^.(9) = i I d'xe''^^^o\T{Mx)4{Omo) , (2) 
where J(x) = J"(x),r/"(x) and J^(x) = J^(x), 

7y"(x) = e''''ql{x)C^5Scix) , 
J«(x) = e-''ulix)Cj,d,{x) , 

= e-''ql{x)Cj,s,{x) , (3) 

the currents J{x) and J^{x) interpolate the scalar and axial- vector diquark states, respec- 
tively, the a, b, c are the color indexes, the C is the charge conjunction matrix, and the 
l^'o) is the nuclear matter ground state. 

In the limit j^o) ~^ |0), we insert a complete set of intermediate "hadronic" states 
with the same quantum numbers as the current operators J{x) and J^{x) into the corre- 
lation functions n(p) and Ii^u[p) to obtain the "hadronic" representation [7J, then isolate 
the ground state contributions from the scalar and axial-vector diquarks, and obtain the 
results: 



= n(,) + + . . . . (4) 

where the pole residues fs and /a are defined as {^\J{^)\S{q)) = fs and {0\J^{0)\A{q)) = 
/a^^i, the is the polarization vector. 

Here we will take a short digression to discuss the application of the QCD sum rules 
in studying the diquark states. In the QCD sum rules, we perform the operator product 
expansion at not so deep Euclidean space, where the approximation of the correlation 
functions by perturbative terms plus some nonperturbative terms makes sense and the 
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contributions from the condensates (or nonperturbative terms) are sizeable. There are 
significant differences between the correlation functions of current operators interpolating 
the diquarks and conventional handrons, we can continue the hadronic correlation func- 
tions to the physical region for the conventional hadrons, but not for the diquarks, as the 
diquarks are non-asymptotic states, there are significant differences between the diquark 
states and conventional hadrons. The one-gluon exchange results in strong attractions in 
the color antitriplet channel 3c, the quark-quark system maybe form quasibound states or 
loosely bound states (diquark states), which are characterized by the correlation length 
L. At the distance / > L, the 3c diquark state combines with one quark or one 3c an- 
tidiquark to form a baryon state or a tetraquark state, while at the distance / < L, the 
3c diquark states dissociate into asymptotic quarks and gluons gradually. We can take 
the diquark state B as an effective colored hadron and the diquark mass as an effective 
quantity, Mjj, ~ , the correlation functions can be continued to the physical region, where 
the quark-quark correlations exist. The transitions two-quarks O diquarks O hadrons are 
not abrupt, the typical correlation lengths L have uncertainties, we have the freedom to 
choose somewhat larger or smaller diquark masses in model-buildings. The correlation 
functions are approximated by a pole term plus a perturbative continuum. 

We use the dispersion relation to express the invariant functions Il{qo, q) in the follow- 
ing form: 



uj-qo 

An(u;, q) = limite_>.o [n('^ + i^^q) — — i^-, q)] ■ (5) 

In the nuclear matter, the corresponding imaginary parts of the spectral densities can be 
expressed as 

An(^, q)=m [F+5{qo - M+) - F_6{qo + M_)] , (6) 

where the Fj- and M± are the pole residues and masses, respectively. In the vacuum limit, 

F± = 'm^^'mX ^^'^ ~ ^S:^A- Thereafter we will use the same notations for the 
masses and pole residues both in the vacuum and in the nuclear matter for simplicity. 

We carry out the operator product expansion in the finite nuclear density at large 
spacelike region q^ <^ 0, and express the invariant functions n(go, q) at the level of quark- 
gluon degrees of freedom as [9l [10] , 

n(go,^) = '^Cn{qo,q){On)pN , (7) 

n 

where the Cniqo^q) are the Wilson coefficients, the in- medium condensates {On)pj^ = 
(^o|0„|^'o) = {O) + pn{0)n at the low nuclear density, the {O) and {0)n denote the 
vacuum condensates and nuclear matter induced condensates, respectively, then take the 
limit = (1,0), = g^, and obtain the imaginary parts of the QCD spectral densities 
according to Eq.(5). One can consult Refs.[9l[T0] for the technical details in the operator 
product expansion. In calculations, we consult the QCD sum rules for the light-flavor 
scalar and axial-vector diquark states in the vacuum, and take the analytical expressions 
of the perturbative terms and the dimension-6 terms from Ref . |il3j . 
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We can match the phenomenological side with the QCD side of the spectral densities, 
and multiply both sides with the weight function ue , then perform the integral j!^^^ dio, 



L 



(8) 



where the loq is the threshold parameter, finally obtain the following two QCD sum rules 
in the nuclear matter: 



r «,(r) 

as{T) 



1 + 



17a, (r) 



e 



T2 



+ o \ _ /Pn + 



TT 



2^2 



rp2 ' 2^2 

167ra, [{qq)l^ + (^g)^^] 

27r2 



(9) 



_m2 



where as{T) 



asiT) 



1 + 



27r 



27r2 [as(M)_ 
8{qUDoq)p^ 8{s''iDos) 

PN 



2ms{qq)pN + 



3mi 
e 

2ms(ss)pjv 



(^1 - e"5*) 



8ms{qiDoiDo)p^ 1 . agGG 56Tras{qq){ss)pj^ 



5r, A = 0.375 GeV, so 



4tt 



9 log 



(10) 

IGeV^. We differentiate Eqs. (9- 



10) with respect to 7^, and obtain two derived QCD sum rules, then eliminate the pole 
residues fs and and obtain the QCD sum rules for the diquark masses. We can replace 
the mass and condensates of the s-quark with the corresponding ones of the g-quark, and 
obtain the QCD sum rules for the ud diquark states. In the limit pN — >■ 0, we obtain the 
corresponding QCD sum rules in the vacuum. The renormalization group improvement 

« 1 at the interval = (0.5 — 1.1) GeV^, while the terms proportional 



factor 



to 



^2 L 



in the derived QCD STim rules have large values, which have significant 



impacts on the masses of the diquark states both m the vacuum and m the nuclear matter. 
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In Fig.l, we plot the values of the -4- 

rp2 



log 



with variations of the Borel parameter 



3 Numerical Results 



In calculations, we have assumed that the linear density approximation is valid at the 
low nuclear density. The input parameters are taken as (g^g)p^ = ^Pn, {s^s)pj^ = 0, 
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Figure 1 : The values of the df 



d 



log 



as(T) "| 9 



with variations of the Borel parameters. 



{m)p. = im) + T^,PN, {ss),, = {ss) + y^^PN, (^)p^ = (^) - (0.65 ± 
0.15)GeV/)7v, {q^iDoq)p^ = O.lSGeVpN, {sUDos)p^ = ^(ss)^^ + 0.02 GeVpjv, + 
ma = 12MeV, = (45 ± 10)MeV, {qHD^iDoq)p^ + ^{q^gsaGq)p^ = 0.031 GeV^TV, 
(sUDoiDos) pj^ + j^{s^gsaGs)pj^ = yO.031 GeV^pTv, {qgsaGq) = ml{qq), (sgsCrGs) = 
ml{ss), (qiDoiDoq) pj^ + \{qgsaGq)p^ = 0.3GeV^/3Ar, {siDoiDQs)p^ + \{sgsaGs)p^ = 
yO.SGeV'^PN, {qgsCrGq)p^ = {qgsaGq)+3.0Ge\''^pN, {sgsCrGs)p^ = {sgsaGs)+y3.0GeV'^pN, 
{q^gs<TGq)p^ = -0.33GeVV7V, {s^gsaGs)p^ = -yO.33 GeV^TV, {qq) = -(0.23 GeV)^, 
{ss) = 0.8{qq), = O.SGeV^, = (0.11 GeV)^, {qq)^^ = f{qq)^^ x (qq),^ + (1 - 
f){qq) X {qq), pn X {^^)pN + (1 ~ X (^•5)5 {qq){ss)pf^ — f{qq)pff X 

(ss)p^ + (1 - f){qq) X (ss), (^^^) = (0.33 GeV)^ / = 0.5 ± 0.5, y = 0.3, = 0.13GeV 
at the energy scale 1^ = 1 GeV [9J . 

In the conventional QCD sum rules [3, there are two criteria (pole dominance and 
convergence of the operator product expansion) for choosing the Borel parameter and 
threshold parameter sq- In this article, we take the pole contributions R as (45 — 80)%, just 
like the ones in our previous works on the heavy, doubly heavy and triply heavy baryon 
states [E], the pole contributions R are defined by 



where the p{s) denotes the QCD spectral densities, the integral over the s can be car- 
ried out analytically, see Eqs.(9-10). In calculations, we observe that larger threshold 
parameters lead to larger Borel windows, and choose the possible smallest Borel windows, 
'^max ~ '^min ~ 0.4 GcV^, to obtain the lowest threshold parameters, and therefore obtain 
the possible lowest masses, which correspond to the largest correlation lengths, the reve- 
lent values are shown explicitly in Table 1. If the same threshold parameters are taken, 
the pole contributions are almost the same in the Borel windows for the QCD sum rules 
both in the vacuum and in the nuclear matter. Thereafter, we will not distinguish the 
pole contributions in the vacuum and in the nuclear matter. On the other hand, the main 
contributions come from the perturbative terms, the operator product expansion are well 



R 




(11) 
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convergent, the two criteria of the QCD sum rules are satisfied. 

Finahy we obtain the numerical values of the masses and pole residues both in the 
vacuum and in the nuclear matter, which are also shown explicitly in Table 1 and Fig. 2. 
The present predictions M^(^(o+) = 0.50 GeV, M^g(Q+) = 0.64 GeV are larger than the 
values Mud{o+) = 0.40 GeV, Mrfg(o+) = 0.46 GeV obtained in Ref.jH], where the deriva- 
tive -TT- does not act on the as(T). From Table 1 and Fig. 2, we can see that including 



the renormalization group improvement factor ^ reduces the diquark masses signif- 

icantly, M^a(o+)-Mud{o+) = 0.14 GeV, Mg,(o+)-Mq^(o+) = 0.13 GeV, M^a(^^+)-M^a{i+) = 
0.11 GeV, Mgs(i+) - Mgs(i+) = 0.12 GeV. Although the factor 



1, the terms 



play an important role in the derived QCD sum 



as(^J.) 
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are rather large compared with the 



associate with the derivative -rr- 

rules. The values of the derivative -A- log 

diquark masses, see Fig.l and Table 1. 

The values of the in- vacuum diquark masses from different theoretical approaches vary 
in a large range, for examples, iV4^(o+) = (0.14 - 0.74) GeV [17], M„rf(o+) = 0.74 GeV, 
Mgs(o+) = 0.88 GeV, M„rf(i+) = 0.95 GeV, M,,(i+) = 1.05 GeV [IS], M^am = 0.82 GeV, 
M,,(o+) = 1.10 GeV, M„rf(i+) = 1.02 GeV, Mg,(i+) = 1.30 GeV [H], M^aio+) = 0.76 GeV, 
^qs{o+) = 0.98 GeV [20j from the Bethe-Salpeter equation with different confining po- 
tentials; M„rf(o+) = 0.710 GeV, Mq^(o+) = 0.948 GeV, M^d{i+) = 0.909 GeV, Mg^(i+) = 
1.069 GeV from a relativistic quark model based on a quasipotential approach in QCD 
[21]; M„rf(o+) = (0.694 ± 0.022) GeV, M,rf(i+) - M„rf(o+) = (0.104 ± 0.042) GeV from 
the lattice QCD [22]; M^d{o+) = (0.42 ± 0.03) GeV, M^a{i+) = (0.94 ± 0.02) GeV [23], 
^ud{o+) = 0.5 GeV [23] from the random instanton liquid model; M„(^(o+) = 0.234 GeV, 
Mud(i+) = 0.824 GeV from the Nambu-Jona-Lasinio Model [25]; M^d(o+) = 0.395 GeV, 
^qs{o+) = 0.590 GeV from the constituent diquark model [26]; etc. One should be careful 
when using them, naively, we expect that they should obey the approximated SUf{3) 
symmetry and the hypersplitting color-spin interaction maybe account for the 0"*" and 1"*" 
diquark mass breaking effects. 

Lattice QCD calculations indicate that the strong attraction in the scalar diquark 
channels favors the formation of good diquarks, the weaker attraction in the axial-vector 
diquark channels maybe form bad diquarks, the energy gap between the light-flavor 
axial- vector and scalar diquarks is about | of the A-nucleon mass splitting, 0.2 GeV 
|27| . which is also expected from the hypersplitting color-spin interaction Tj • Tjcii ■ aj 
[H [5]. The coupled rainbow Dyson-Schwinger equation and ladder Bethe-Salpeter equa- 
tion also indicate such an energy hierarchy [TP]. In the present work, the central values 
have the energy gaps M„rf(i+) - M^a(o+) = 0.17GeV, Mg^(i+) - Mg3(o+) = 0.15 GeV, 

^ud{i+) ~ ^ud{o+) = 0.20 GeV, Mg^(x+) ~ ^qs{o+) — 0.16 GeV, which are consistent with 
predictions of the lattice QCD and Bethe-Salpeter equation. If we neglect the uncertain- 
ties, the SUf{3) breaking effects for the masses of the scalar and axial-vector diquark states 
are Mgs{o+) - M„d(o+) = 0.13 GeV, Mg^(i+) - M^dii+) = 0.11 GeV, Mgs(o+) - A4^(o+) = 
0.14 GeV, Mg^(x+) ~ ^ud(i+) = 0.10 GeV, respectively, which are consistent with the naive 
expectation — mg 0.13 GeV. 
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pole 


M(GcV) 


/(GeV^) 


uuyu ) 


n c; n Q /I 9 
u.o — u.y/ i.z 


(_40 — to) /O 


U.D4i It U.UD [U.Do It U.UDJ 


It U.U± I It U.UiDJ 


qs{0+) 


0.6-1.0/1.4 


(47 - 74)% 


0.77 ± 0.04 [0.79 ± 0.04] 


0.313 ± 0.013 [0.312 ± 0.012] 


ud{l+) 


0.6 - 1.0/1.6 


(48 - 76)% 


0.81 ±0.06 [0.87 ±0.05] 


0.228 ±0.016 [0.233 ±0.012] 


qs{l+) 


0.7-1.1/1.8 


(49 - 74)% 


0.92 ± 0.04 [0.95 ± 0.03] 


0.269 ±0.011 [0.269 ±0.010] 


ud{0+) 


0.5 - 0.9/1.2 


(45 - 77)% 


0.50 ± 0.05 [0.54 ± 0.04] 


0.246 ± 0.010 [0.243 ± 0.009] 


qs{0+) 


0.6 - 1.0/1.4 


(47 - 74)% 


0.64 ± 0.03 [0.65 ± 0.03] 


0.286 ±0.007 [0.283 ± 0.007] 


ud{l+) 


0.6 - 1.0/1.6 


(48 - 77)% 


0.70 ± 0.04 [0.74 ± 0.04] 


0.209 ±0.010 [0.210 ±0.008] 


qsil+) 


0.7-1.1/1.8 


(49 - 75)% 


0.80 ± 0.03 [0.83 ± 0.02] 


0.244 ± 0.006 [0.242 ± 0.006] 



Table 1: The Borel parameters, threshold parameters, pole contributions, masses and 
pole residues of the light-flavor diquark states, the wide-hat ^ denotes the renormalization 



group improvement factor 
nuclear matter. 



is included, and the bracket denotes the values in the 



Prom Table 1, we can see that the diquark states in the nuclear matter have larger 
masses than the corresponding ones in the vacuum, the in-medium effects lead to the 
mass-shifts (5M„^(o+) = 0.04 GeV, (5Mg^(o+) = 0.02 GeV, 5M„d(o+) = 0.04 GcV, 5Mqs{o+) = 
0.01 GeV, (5M„d(i+) = 0.06 GeV, 5Mg^(i+) = 0.03 GeV, dM^d{i+) = 0.04 GeV, 6Mqs{i+) = 
0.03 GeV. Although the diquark masses have uncertainties originate from the Borel param- 
eters, the mass-shifts survive approximately as the uncertainties are canceled out with each 
other, see Fig. 2, the quark-quark correlation lengths are reduced slightly in the nuclear 
matter. Compared with the ud diquark states, the qs diquark states have much smaller 
mass-shifts, which attributes to the condensates of the s-quark obtain much smaller mod- 
ifications than the corresponding ones of the u and d-quarks. On the other hand, we can 
see that the pole residues in the vacuum and in the nuclear matter approximately have 
the same values, as the dominating contributions come from the perturbative terms. 



4 Conclusion 

In this article, we study the light-flavor scalar and axial-vector diquark states in the vac- 
uum and in the nuclear matter using the QCD sum rules in an systematic way. The 
predicted diquark masses in the vacuum obey the flavor SUf{S) symmetry approximately, 
and the 0"*~ and 1"*" diquark mass breaking effects are consistent with the lattice calcula- 
tions. The diquark states in the nuclear matter have larger masses than the corresponding 
ones in the vacuum, the quark-quark correlation lengths are reduced slightly in the nu- 
clear matter. We can take the diquark masses as basic parameters and perform many 
phenomenological analysis. 
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Figure 2: The masses M of the hght-flavor diquark states with variations of the Borel 
parameters T^, the A, B, C and D correspond to the ud{0~^), qs{0~^), ud{l^) and gs(l+) 
diquark states respectively. The a, A and /3, r denote the values from the QCD sum 
rules in the nuclear matter and in the vacuum respectively; while the A, r and a, (3 de- 



note the renormalization group improvement factor 
respectively. 
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9 . 



is included and not included 
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